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Abstract 

This pearl aims to demonstrate the ideas of wholemeal and projec- 
tive programming using the Towers of Hanoi puzzle as a running 
example. The puzzle has its own beauty, which we hope to expose 
along the way. 

Categories and Subject Descriptors D.1.1 [Programming Tech- 
niques]: Applicative (Functional) Programming; D.3.2 [Program- 
ming Languages]: Language Classifications — Applicative (func- 
tional) languages, Haskell; D.3.3 [Programming Languages]: 
Language Constructs and Features — Frameworks, patterns, recur- 
sion 

General Terms Algorithms, Design, Languages 

Keywords Towers of Hanoi, wholemeal programming, projective 
programming, Hanoi graph, Sierpinski graph, Sierpihski gasket 
graph, Gray code 

1. Introduction 

Functional languages excel at wholemeal programming, a term 
coined by Geraint Jones. Wholemeal programming means to think 
big: work with an entire list, rather than a sequence of elements; de- 
velop a solution space, rather than an individual solution; imagine 
a graph, rather than a single path. The wholemeal approach often 
offers new insights or provides new perspectives on a given prob- 
lem. It is nicely complemented by the idea of projective program- 
ming: first solve a more general problem, then extract the interest- 
ing bits and pieces by transforming the general program into more 
specialised ones. This pearl aims to demonstrate the techniques us- 
ing the popular Towers of Hanoi puzzle as a running example. This 
puzzle has its own beauty, which we hope to expose along the way. 



if we arrange them in a circle. Initially, the discs are placed on one 
peg in decreasing order of diameter. The task is then to move the 
disks, one at a time, to another peg subject to the rule that a larger 
disk must not be placed on a smaller one. 

This restriction implies that a configuration can be represented 
by a list of pegs: the first element determines the position of the 
largest disc, the second element the position of the second largest 
disc, and so forth. Consequently, there are 3" possible configura- 
tions where n is the total number of discs. Lucas' original puzzle 
contained 8 discs. The instructions of the puzzle refer to an old 
Indian legend, attributed to the French mathematician De Parville, 
according to which monks were given the task of moving a total 
of 64 discs. Since the day of the world's creation, they transfer the 
discs, one per day. According to the legend, once they complete 
their sacred task, the world will come to an end. 

Now, taking a wholemeal approach, let us first develop the 
big picture. The set of all configurations together with the set 
of legal moves defines a graph, which turns out to enjoy a nice 
inductive definition. If there are no discs to move around, the graph 
consists of a singleton node: o. For the inductive step we reason as 
follows: the largest disc can only be moved, if all the smaller discs 
reside on one other peg. The smaller discs, however, can be moved 
independent of the largest one. As the largest disk may rest on one 
of three pegs, the graph forn -I- 1 discs consequenfly incorporates 
three copies of the graph for n discs linked by three edges. The 
diagram in Fig. 1 illustrates the construction. The graph has the 
shape of a triangle; the dashed lines indicate the sub-graphs (for 
n = 0 the sub-graphs collapse to singleton nodes); the three solid 
lines connect the sub-graphs. The inductive construction shows that 
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2. The Hanoi graph 

The Towers of Hanoi puzzle was invented by the French number 
theorist Edouard Lucas more than a hundred years ago. It consists 
of three vertical pegs, on which discs of mutually different diame- 
ters can be stacked. For reference, we call the pegs A, B and C and 
let a, b and c range over pegs. 

data Pes = A I B I C 

1 own a version of the puzzle where the pegs are arranged in a row. 
However, the mathematical structure of the puzzle becomes clearer. 
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Figure 1. Inductive construction of the Hanoi graph. 

the graph is planar: it can be drawn so that no edges intersect. Fig. 2 
displays the graph for 4 discs. To reduce clutter, the peg of the 
largest disc is always written in the centre of the respective sub- 
graph, with the size of the font indicating the size of the disc. Can 
you find the configuration [B, A, B, c]? The comers of the triangle 
correspond to perfect configurations: all the discs reside on one 
peg. The example graph shows that every configuration permits 
three different moves, except for the three perfect configurations, 
where only two moves are possible. 



Figure 2. The Hanoi graph for 4 discs. 



Let us turn our attention to the layout of the sub-graphs. The fol- 
lowing notation proves to be useful: the arrangement z denotes 
a permutation of x, y and z, which by assumption are pairwise dif- 
ferent. Using a'c to indicate the position of the largest disc — the 
largest disc in the left triangle resides on a and so forth — we ob- 
serve that if the comers of the graph are a ' c , then the comers of the 
sub-graphs are a°fr, c^'a and t^c, respectively. Using this observa- 
tion, we can capture the informal description of the constmction as 
a pseudo-Haskell program. 

graphg {J c) = o 

graph^^-^ („'c) = b< graph,^ {J' a) 

I \ 
a < graph^^ {a" t) — c < graph^ 

The function graph^ maps an arrangement a'c to an undirected 
graph, whose vertices are lists of pegs of length n. The notation 
a < g means prepend a to all the vertices of g; o denotes a 
singleton node labelled with the empty list. We leave the type of 
graphs unspecified. If the type were a functor, then a < g would be 
fmap {a :) g. The call graph^ (a'^c;) yields the graph in Fig. 2. 

A few observations are in order. The definition of graph^ im- 
plies that the graph has 3" nodes (ao = 1, a„+i = 3 • a„) and 
_ 3~) ^2 edges (ao = 0, On+i = 3 ■ a„ + 3). Furthermore, the 
length of a side of the triangle is 2" — 1 (ao — 0, a„+i = 2- a„ -|- 1). 
Since there are only 3! = 6 permutations of three different items, 
the graph contains at most six different mini-triangles: a^c, c^b, 
b'^a, a'^b, c^a and e^c- Note that the first three arrange the 
pegs clockwise and the last three coimterclockwise. Inspecting the 
definition of graph ^, we see that the direction changes with every 
recursive step. The diagram in Fig. 3 illustrates the change of di- 
rection. This observation implies that the recursion pattern is quite 
regular: at any depth there arc only three different recursive calls. 
Fig. 4 visualises the call stracture using three different colours. 




Figure 3. Change of direction. 



3. Towers of Hanoi, recursively 

Solving the Towers of Hanoi puzzle amounts to finding a shortest 
path in the corresponding graph. Clearly, the shortest path between 
two comers is along the side of the triangle. Projecting onto the 
lower side, we transform graph^ to 

hanoi'o {a'' c) = [[]] 

hanoi'n+i (a'c) = a < hanoi'^ (a^b) + c < hanoi'^ {b°'c) ■ 

Now, a < X is shorthand for /map (a :) x. The call hanoi'„ (a^'c) 
returns a list of configurations solving the problem of 'moving n 
discs from a to c using b\ 

Instead of configurations (vertices of the graph) we can alterna- 
tively return a list of moves (corresponding to edges). 

hano'H) {Jc) = [] 

hanoin+i (o''c) = hanoi„ [a^b) + [(a, c)] + hanoi„ (b^c) 



Figure 4. Call structure of hanoi at recursion depth 4. 



The pair (a, c) represents the instruction 'move the topmost disc 
from a to c'. Here is the sequence of moves for n = 4: 

:> hanoii (a^c) 

[(A, B), (A, C), (B, C), (A, B), (C, A), (C, B), (A, B), (A, C), 
(B, C), (B, A), (C, A), (B, C), (A, B), (A, C), (B, C)] . 

Note that the lower part of the list can be obtained from the upper 
part via a clockwise rotation of the pegs: a'^h becomes (,"c- 

There is at least one further variation: instead of an arrangement 
one can pass a source and a target peg. 

hanoio a c — [] 
hanoin+i o, c 

— hanotn a (a ± c) + [(a, c)] + hanoin {a ± c) c 

The function ±, which determines 'the other peg', is given by 



A±A = A 
B±A = C 

C± A = B: 



A. 
B . 

C . 



B = C 
B = B 
B = A 



A. 
B. 
C. 



C = B 
C = A 
C = C 



We will find some use for ± later on. For the moment, we just 
note that the operation is commutative and idempotent, but not 
associative. 

4. Towers of Hanoi, parallelly 

Imagine that the monastery always accommodates as many monks 
as there are discs. The tallest monk is responsible for moving 
the largest disc, the second tallest monk for moving the second 
largest disc, and so forth. Can we set up a work schedule for the 
monastery? 

Inspecting Fig. 2, we notice that, somewhat unfairly, the small- 
est monk is the busiest. Since the smallest triangles correspond to 
moves of the smallest disc, he is active every other day. We can ex- 



tract his work plan from hanoin by omitting all the moves, except 
when n equals 1. 

cyckg (Jc) = [] 

cyde-^ (Jc) = [{a,c)] 

cyde„+i (a'c) = cj/c/e„ + cycle^ 

The function is called cycle, because the smallest disc cycles 
around the pegs: in the recursive call it is moved from a to 6 
and then from b to c. Whether it cycles clockwise or counterclock- 
wise depends on the parity of n — the direction changes with every 
recursive call of cycle. 

Of course, the smallest disc is by no means special: all the discs 
cycle around the pegs, albeit at a slower pace and in alternating 
directions. In fact, hanoi satisfies the following 'fractal' property: 

hanoin+i (o''c) = cycle^j^-^ (J c) Y hanoin (a'c) , (1) 

where Y denotes the interleaving of two lists. 

[] Y bs = bs 

(a : as) Y bs = a : {bs Y as) 

The fractal property suggests the following alternative definition of 
hanoi, which has a strong parallel flavour. 



phanoig (a'c) 
phanoi^^-^ (Jc) 



cycle^^j^ (a'c) Y phanoi^ (Jc) 



In words, the smallest monk starts to work on the first day; he is 
active every second day and moves his disc, say, clockwise around 
the pegs. The second smallest monk starts on the second day; he 
is active every fourth day and moves his disc counterclockwise. 
And so forth. Actually, only the smallest monk must remember the 
direction; for the larger discs there is no choice, as one of the other 
two pegs is blocked by the smallest disc. 



There is an intriguing cross-connection to binary numbers: the 
activity diagram of the monks (Which monk has to work on a given 
day?) 

discso = [] 

discsn+i = discsn + [n] + discs„ 
yields the binary carry sequence or ruler function as n goes to 
infinity (Hinze, 2008). This sequence gives the number of trailing 
zeros in the binary representations of the positive numbers, most 
significant bit first. Or put differently, it specifies the running time 
of the binary increment. We will make use of this observation later 
on. 

The fractal property (1) enjoys a simple inductive proof, which 
makes essential use of the following abide law. If xi and X2 are of 
the same length, then 

[xi + yi) Y (3>2 + 2/2) = [xi Y X2) + (2/1 Y 2/2) . 

The basis of the induction is straightforward. Here is the inductive 
step. 

cyc/e„_,_2 [a'' c) Y hanoi„+i (Jc) 
= { definition of cycle and definition of hanoi } 

{cyclen+i U^h) + cycle 

n + l {i <:)) 

Y {hanoin (a%) + [{a, c)] + hanoi„ {b" c)) 
= { abide law } 

(c2/cfe„_|_i {a^h) Y hanoin (0%)) + [(a, c)] 
+(c?/c;e„+i Y hanoin ((."c)) 
= {ex hypothesi } 

hanoin+i (a^t) + [{a, c)] + hanoin+i (b^c) 
= { definition of hanoi } 

hanoin+2 {J c) 

The abide law is applicable in the second step, because the two 
lists, cycle^j^-i {a'' b) and hanoin (a'^b) + [(a, c)], have the 
same length, namely, 2". Speaking of the length of lists, note 
that 2"+^ — 1 = Z)"=o 2' ^ simple consequence of the fractal 
property. 

5. When will the world come to an end? 

Many visitors come to the monastery. Looking at the configuration 
of discs, they often wonder how many days have passed since the 
creation of the world. Or, when will the world come to an end? 

We can answer these questions by locating the current config- 
uration in the Hanoi graph. If we use as positions the three digits 
0^1 — that is, 0 for the left triangle, 2 for the upper triangle and 1 
for the right triangle — then we obtain the answer to the first ques- 
tion in binary. Fig. 5 displays the Hanoi graph for 4 discs suit- 
ably re-labelled — this graph is also known as the Sierpinski graph. 
The 2'^ nodes on the lower side of the triangle, and only those, are 
labelled with binary numbers. Consequently, if the current position 
contains a 2, we know that the monks have lost track. (To answer 
the second question, we use as positions 1^0 instead of 0^1. If we 
are only interested in the distance to the final configuration, we 
simply replace the digit 2 by a 1.) 

POS (a^) [] = [] 
POS (a'c) (p ■■ ps) 

I p o = 0 : pos (a^b) ps 
\ p == b =2: pos {c'' a) ps 
\ p --^ c =1 : pos (t^c) ps 

For instance, pos (a^c) [A,B,C,A] yields [0,1,0,1], the bi- 
nary representation of 5, most significant bit first. The function 
pos (a'c) defines a bijection between {A, B, C}" and {0, 2, 1}" 



for any given initial arrangement a c • This arrangement can be seen 
as representing the bijective mapping {a 1-^ 0, b i—> 2, c 1-^ 1}. Al- 
ternatively, we can use an arrangement, say, ;*r as a representation 
of the 'inverse' mapping {A i—^ l,Ti i—> t,C l—^ r} obtaining the 
following slightly more succinct variant of pos. 

pos' UK) [] = [] 

pos' {A: ps) = I : pos' (^t) ps 

pos' (I'r) {B : ps) = t : pos' (r'i) ps 
pos' {i\) (C : ps) = r : pos' (tV) ps 

The two variants are related by pos (a^c) = pos' (0^1). The 
latter definition is particularly easy to invert. 

conf (Jc) [] = [] 

conf (aK) {0 : ps) = a : conf (a^b) ps 
conf {a''c) {2 : ps) = b : conf (c'a) ps 
conf (Jc) (1 : ps) = c : conf (b^c) ps 

The call conf (a^c) [0, 1, 0, 1] yields [A, B, C, A], the configu- 
ration we obtain after 5 days. The functions pos' and conf are actu- 
ally identical, if we identify A with 0, B with 2 and C with 1. Then 
pos' is an involutive graph isomorphism between Hanoi graphs and 
Sierpinski graphs of the same order. 

If the monks have lost track — a 2 appears in the answer list — 
then we can use the idea of locating a configuration in the Hanoi 
graph to determine the shortest path to the final configuration c". 
This leads to the following generalised version of hanoi. 

ghanoig [] {a\) ^ [] 

ghanoin+i (p : ps) (Jc) 

I p a = ghanoin ps (a^b) + [{a, c)] + hanoin {b" c) 
I p == 6 = ghanoi^ ps (b^a) + [{b, c)] + hanoin {a*" c) 
I p c = ghanoin P^ ("''c) 

(The index n is actually redundant: it always equals the length of 
the peg list.) Depending on the location of p, we either walk within 
the left triangle and then along the lower side, or within the upper 
triangle and then along the right side, or within the right triangle. 
The reader should convince herself that ghanoi^ ps (a'c) indeed 
yields the shortest path between ps and c". Note that the first two 
equations are perfectly symmetric, in fact, ghanoi^ ps {J" c) = 
ghanoin P^ (b'^c)- As an aside, if we fuse ghanoi with length, 
then we obtain a function that yields the distance to the final 
configuration. 

6. Towers of Hanoi, iteratively 

The generalised version of the puzzle — moving from an arbitrary 
configuration to a perfect configuration — serves as an excellent 
starting point for the derivation of an iterative solution, that is, a 
function that maps a configuration to the next move or to the next 

configuration. 

The next move is easy to determine: we fuse ghanoi with the 
natural transformation 

first :: [a] —> Maybe a 

first [] — Nothing 

first {a : as) = Just a 

that maps a list to its first element. We obtain 

moveo [] (a'c) = Nothing 

move„+i (p : ps) (a'c) 

\ p -- a = moven ps (a^b) + Just{a, c) + first{hanoin(b°' c)) 
I p == 6 = moven ps (b^a) + Just{b, c) + first {hanoin {a'' c)) 

I p C = mOVCn ps {a''c) ■ 



Figure 5. The Sierpinski graph of order 4. 



The operator + is overloaded to also denote concatenation of 
Maybe values (+ is really the rnplus method of MonadPlus). 

(+) :: Maybe a — > Maybe a — > Maybe a 

Nothing + m = m 

Just a + m — Just a 

We can simplify the definition of move drastically: since Just a + 
m = Just a, the calls to hanoi can be eUminated; because of 
that the index n is no longer needed. Furthermore, the first two 
equations can be unified through the use of ±, the operator that 
determines 'the other peg'. Finally, the second argument, a'' c, can 
be simplified to c, the target peg. 

move [] c = Nothing 
move {p : ps) c 

\ p ^ c = move ps {p ± c) + Just {p, c) 

\ p == c = move ps c 

Here is a short interactive session that illustrates the use of move. 

> move [A, B, B, A, C, C, C, B] C 
Just (B, C) 

> move [A, B, B, A, C, C, C, C] C 
Just (A, B) 

> move [C, C, C, C, C, C, C, C] C 
Nothing 

Since the pair (B, C) means 'move the topmost disc from peg B 
to peg C, the configuration following [A, B, B, A, C, C, C, B] is 
[A, B, B, A, C, C, C, C]. Incidentally, if we start with the initial 
peg list [A, A, A, A, A, A, A, A] and target C, we obtain these 
configurations after 110 and 111 steps. 

The function m.ove implements a two-way algorithm: on the 
way down the recursion it calculates the target peg for each disc 
(using ±); on the way up the recursion it determines the smallest 
disc that is not yet in place (using +). The following table makes 



the target pegs explicit. The rows labelled c, lists the target peg for 
each disc: the first, user- supplied target is co = C, the next target 

is ci = po -L Co, and so forth. 

i 0 12 :j 1 5 () 7 

110 Pi A B B A C C C B~ 

a CBBBCCCC A 

111 Pi A B B A C C C C 

c, CBBBCCCCC 

The smallest disc is not in place, so it is moved from B to C. In the 
next round, disc 3 has to be moved from A to B. 

The smaller discs are moved more frequently, so it is actually 
prudent to reverse the list of pegs so that the peg on which the 
smallest disc is located comes first. The situation is similar to the 
binary increment: visiting the digits from least to most significant 
bit is more efficient than the other way round. Let us assume for 
the moment that we know the 'last target', the value of c that is 
discarded in the first equation of move (the pegs that stick out in 
the example above). Since ± is reversible, p± c = c' iff c = p± c', 
we can reconstruct the previous target c from the next target c'. 
The following variant of move makes use of this fact — the suffix 
'i' indicates that the input list is now arranged in increasing order 
of diameter. 

movei [] c' = Nothing 

movei {p : ps) c' 

\ p {p ± c') = Just (p, p ± c') + movei ps {p ± c') 
\ p -- {p ± c') = movei ps c' 

We consistently changed the second argument of move to reflect 
the fact that we compute the previous from the next target and 

additionally replaced the remaining occurrences of c by p ± c'. 
Again, we can simplify the code: p ^ (p ± c') is the same as 
p ^ c' . Applying Just a + m = Just a once more, we can 
eliminate the first recursive call to movei so that the function 



stops as soon as the smallest displaced disc is found — this was the 
purpose of the whole exercise. We obtain 

movei [] c' = Nothing 

movei (p : ps) c! 

\ p ^ c' = Just {p,p J- c') 

\ p -- c' = movei ps c' . 

In a nutshell, movei determines the smallest disc that does not 
reside on the last target. 

So, for an iterative version of hanoi we have to maintain two 
pieces of information: the current configuration and the current 'last 
target'. It remains to determine the initial last target and how the 
last target changes after each move. If the Ust of pegs is given in the 
original decreasing order, then we can transform move to 

last [] c = c 

last (p : ps) c = last ps {p ± c) , 

which yields the last target. It is not hard to see that last is an 
instance of the famous foldl: last ps c = foldl (±) c ps. If 
we reverse the list, we simply have to replace foldl by foldr addi- 
tionally using the fact that ± is commutative: foldr (±) c ps = 
foldl (±) c {reverse ps). 

Next, we augment movei so that it returns the next configura- 
tion instead of the next move and additionally the next 'last target'. 

step ([], c') = Nothing 

step (p : ps, c') 

\ p c' = Just {{p ± c') : ps, p ± c') 
\ p == c' = do { {ps' , c) <— step {ps, c'); 

return {p : ps' , p ± c)} 

Consider the second equation: p is moved to p± c'. After the move, 
the disc resides on the target peg. Consequently, the next target is 
also p ± c' — ^recall that ± is idempotent. This target is then updated 
in the third equation 'on the way back' mimicking move's mode of 
operation. The function step runs in constant amortised time, since 
it performs the same number of steps as the binary increment — 
recall that the activity diagram of the monks coincides with the 
binary carry sequence. 

The next last target can, in fact, be easily calculated by hand. 
Consider the following two subsequent moves (as usual, a, b and c 
are pairwise different). 

n n — 1 n — 2 10 

• • • a c c ■ ■ ■ c c 

■ ■ ■ h c c •■•ccc 

• • • 6 c c • • • c c 
•••6 h a ■ ■ ■ h a b 

Assume that the first configuration ends with an even number of cs. 
The topmost disc of a is then moved to b. The new succession 
of target pegs consequently alternates between b and a: since the 
number of cs is even, the new last target is b\ for an odd number, 
it is a. So, the monks can be instructed as follows: determine the 
smallest disc that is not on c. Transfer it from a to b. If the disc's 
number is even, the new last target is 6; otherwise, it is a. 

If we solve the original puzzle, that is, if the configurations he 
on one of the sides of the triangle, then the next last target is even 
easier to determine: like the discs, it cycles around the pegs. If the 
pegs are arranged a^c and we move the discs from A to C, then 
the last target moves counterclockwise for an even number of discs 
and clockwise for an odd number. 



Summing up, we obtain the following iterative implementation 
for solving the generalised Hanoi puzzle. 

ihanoi ps c = map fst {iterate step {ps, foldr (±) c ps)) 

iterate :: {a — > Maybe a) — > (a ^ [a]) 
iterate fx = x : case / a; of { Nothing —> []; 

Just x' — > iterate fx'} 

7. Longest paths and Sierpinski's triangle 

So far we have considered shortest paths in the Hanoi graph. Since 
the destruction of the world hangs in the balance, as a gift to 
future generations, we might want to look for the longest path. 
In the following variant of hanoi the largest disc makes a detour. 
(According to the Indian legend, the temple is actually in Bramah 
rather than in Hanoi, hence the name of the function.) 

bramaho {a'' c) = [] 

bramah„+i {a^c) = bramahn (a'c) + [{a, b)] 
+ bramahn (c'a) + [{b, c)\ 
+ bramahn (o'c) 

The largest disc is first moved from a to b, and then from b to c. 
Since bramahn returns 3" — 1 moves (ao = 0, an+i = 3 ■ an + 2), 
we have actually found a Hamiltonian path. The path has another 
interesting property: if the pegs are arranged in a row, a b c, then 
discs are only moved between adjacent pegs. 

The Hamiltonian path for four discs is displayed in Fig. 6. The 
picture is quite appealing. Actually, if we move the sub-triangles 
closer to each other so that the corners touch, we obtain a nice frac- 
tal image. Fig. 7 shows the result for 7 discs. The corresponding 
graph is known as the discrete Sierpihski gasket graph.' The pic- 
ture has been drawn using Functional Metapost's turtie graphics 
(Korittky, 1998). 

curveo d = forward & turn {2 * d) & forward 
curve„+i d = turn d & curven {—d) & turn d & curve„ d 
& turn d & curven {—d) & turn d 

The command forward moves the turtle one step forward, turn d 
turns the turtle clockwise by d degrees, and & sequences two 
turtle actions. Since turtle graphics is state-based — the turtle has 
a position and faces a direction — ^recursive definitions typically 
maintain an invariant. To draw the 'triangle' o'c, we start at a 
looking at b and stop at c looking away from b. The overall change 
of direction is twice the second argument of curve, which for 
equilateral triangles is either 60° or —60°. 

The curve is closely related to Sierpinski's arrowhead curve. In 
fact, both curves yield Sierpinski's triangle as n goes to infinity. 
As an aside, Sierpinski's triangle is a so-called /racfa/ curve: it has 
the Hausdorff dimension log 3/ log 2 = 1.58496, as it consists of 
three self-similar pieces with magnification factor 2. 

8. Gray code 

The function bramah enumerates the configurations {A, B, C}" 
changing only one peg at a time. In other words, the succession of 
configurations corresponds to a ternary Gray code! To investigate 
the correspondence a bit further, here is a version of bramah that 
returns a Ust of configurations, rather than a Ust of moves. 

bramah'o (o'c) = [[]] 
bramah'n+i (a'c) = a < bramah'n (a'c) 
+ 6 < bramah'n (c''a) 
+ c < bramah'n (a'c) 



The term gasket graph is actually not used consistently. Some authors 
call the counterpart of the Hanoi graph the Sierpidski gasket graph, and 
its contracted variant the Sierpinski graph. 



Figure 6. A Hamiltonian path in tiie Hanoi graph for 4 discs. 




Figure 7. A Hamiltonian path in the Hanoi graph for 7 discs. 



There are two standard ternary Gray codes: the so-called modular 
code, the digits vary 012 1 201 1 120 • • •, and the reflected code, the 
digits vary 012|210|012---. The definition above yields the latter, 
as the discs are only moved between adjacent pegs. In fact, we have 
bramah'n {c'' a) = reverse (bramah'n (a'' c))- Using this property 
on the second recursive call, we can simplify bramah'„ (0^2) 
somewhat. 

gray 3 =[[]] 

gray3^_^i = 0 < gray3^ + 1 < reverse gray3^ + 2 < gray 3 „ 
For comparison, here is the definition of the Gray binary sequence. 
gray2o =[[]] 

gray2^^i = 0 < gray2^ + 1 < reverse gray2^ 

Actually, the binary Gray code is also hidden in the Tower of 
Hanoi puzzle. We only have to work with a different configuration 
space: instead of {A, B, C}" we use {0, 1}" keeping track whether 
a disc has been moved an even or an odd number of times. The 
initial configuration A A" becomes 00", the final configuration 
CC" becomes 10". Since the ith disc is moved 2' times, all the 
discs make an even number of moves, except for the largest, which 
makes a single move. 

We can easily adapt ihanoi to generate binary Gray code. We 
take as a starting point the first definition of movei, slightly mod- 
ified so that the next configuration is returned instead of the next 
move (we also applied the Just a + m = Just a simplification). 

configi [] c — Nothing 

configi (p : ps) c 

I p 7^ (p ± c') = Just {p ± c' : ps) 

\ p -- {p ± c') = p < configi ps (p ± c') 

Now, the Gray code equivalent of ± is the Boolean operation 
exclusive or, that is, inequality of Booleans. This implies that the 
last target c' corresponds to a parity bit. Thus, configi becomes 
(the code uses Booleans rather than bits) 

codei [] P = Nothing 

codei {b : bs) p 

\ b -/- {b -/- p) = Just ((6 -/- p) : bs) 
\ b ~ {b ^ p) = b < codei bs {b ^ p) . 

Again, we can simplify the code a bit: inequality and equality of 
Boolean values are associative (Backhouse and Fokkinga, 2001), 
so 6 7^ (6 p) is simply p. Using False ^ b = b and 
True 6 = -1 6, we obtain 

codei [] p = Nothing 

codei {b : bs) p 

I p = Just (-1 b : bs) 

\ -I p = b <3 codei bs b . 

In words: we traverse the list p : 6s up to the first 1; the following 
bit, if any, is flipped. 

As for ihanoi, we have to maintain two pieces of information: 
the current Gray code and the current parity bit. The latter is easy 
to update: it is flipped in each step. Sunmiing up, we obtain the 
following Gray code generator. 

igray bs = map fst {iterate step {bs,foldr {^) True bs)) 
where step {bs,p) = do { bs' ^ codei bs p; 

return (6s', -■ p) } 

This is, in fact, the functional version of Knuth's Algorithm G 
(2005). 

9. Conclusion and further reading 

We have come to the end of the Tour D' Hanoi. In the spirit of the 
wholemeal approach we started with an inductive definition of the 



Hanoi graph. From that we derived a series of programs evolving 
around the Tower of Hanoi theme. Knowing the big picture was 
jolly useful: for instance, calculating the number of moves could 
be reduced to the problem of locating a configuration in the graph. 
Projective program transformations are abundant: hanoi is derived 
from graph by projecting onto the lower side of the graph, cycle is 
derived from hanoi by mashing out the moves of the larger discs, 
and so forth. 

The transformations could be made rigorous within the Algebra 
of Programming framework (Bird and de Moor, 1997). Occasion- 
ally, this comes at an additional cost. For instance, to derive cycle 
from hanoi we would additionally need the disc number, which is 

not present in hanoi 's output. 

A lot is left to explore. There are literally hundreds of papers 
on the subject: Paul Stockmeyer's comprehensive bibliography has 
a total of 369 entries (2005). From that bibliography I learned that 
my definition of the Hanoi graph is not original: Er introduced it to 
analyse the complexity of the generalised Tower of Hanoi problem 
(1983). The original instructions of the game already alluded to the 
recursive procedure for solving the puzzle. It has been used since 
to illustrate the concept of recursion. The parallel version — usually 
classified as iterative — is due to Buneman and Levy (1980). Back- 
house and Fokkinga (2001) show that each disc cycles around the 
pegs exploiting the associativity of equivalence. To the best of the 
author's knowledge the iterative, or stepwise variant is original. The 
connection to Gray codes has first been noticed by Gardner (1972). 
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